Adiabatic Domain Wall Motion and Landau-Lifshitz Damping 
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Recent theory and measurements of the velocity of current-driven domain walls in magnetic 
nanowires have re-opened the unresolved question of whether Landau-Lifshitz damping or Gilbert 
damping provides the more natural description of dissipative magnetization dynamics. In this 
paper, we argue that (as in the past) experiment cannot distinguish the two, but that Landau- 
Lifshitz damping nevertheless provides the most physically sensible interpretation of the equation 
of motion. From this perspective, (i) adiabatic spin-transfer torque dominates the dynamics with 
small corrections from non-adiabatic effects; (ii) the damping always decreases the magnetic free 
energy, and (iii) microscopic calculations of damping become consistent with general statistical and 
thermodynamic considerations. 
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I. BACKGROUND 

Experiments designed to study the effect of elec- 
tric current on domain wall motion in magnetic 
nanowires show that domain walls move over large 
distances with a velocity proportional to the applied 
current ±&2d£& L 8 1 9 1 1 Most theories ascribe this behav- 
ior to the interplay between spin-transfer (the quantum 
mechanical transfer of spin angular momentum between 
conduction electrons and the sample magnetization) and 
magnetization damping of the Gilbert typeJ^ Contrary 
to the second point, we argue in this paper that Landau- 
Lifshitz damping^ provides the most natural descrip- 
tion of the dynamics. This conclusion is based on the 
premises that damping should always reduce magnetic 
free energy and that microscopic calculations must be 
consistent with statistical and thermodynamic consider- 
ations. 

Theoretical studies of current-induced domain wall 
motion typically focus on one-dimensional models where 
current flows in the x-direction through a magnetization 
M(a;) = MM(l). When M is constant, the equation of 
motion is 



M 



-7M x H + N ST + D. 



(1) 



The precession torque — 7M x H depends on the gyro- 
magnetic ratio 7 and an effective field /xqH = — 6F/6M. 
which accounts for external fields, anisotropies, and 
any other effects that can be modelled by a free en- 
ergy F[M] (fio is the magnetic constant). The spin- 
transfer torque Nst is not derivable from a poten- 
tial, but its form is fixed by symmetry arguments 
and model calculations^^^^Li 8 .^^^ A local 
approximation^ (for current in the x-direction) is 



The first term in ^ occurs when the spin current follows 
the domain wall magnetization adiabatically, i.e., when 
the electron spins remain largely aligned (or antialigned) 
with the magnetization as they propagate through the 
wall. The constant v is a velocity. If P is the spin polar- 
ization of the current, j is the current density, and /is is 
the Bohr magneton, 



eM ' 



(3) 



The second term in © arises from non-adiabatic effects. 
The constant is model dependent. 

The damping torque D in |T]) accounts for dissipa- 
tive processes, see [23| for a review. Two phenomenolog- 
ical forms for D are employed commonly: the Landau- 
Lifshitz formic with damping constant A, 

D L = —AM x (M x H) , (4) 

and the Gilbert formii with damping constant a, 



D fi = aM x M. 



(5) 



N 



ST 



d x M - /3M x d x M 



(2) 



The difference between the two is usually very small 
and almost all theoretical and simulation studies of 
current-induced domain wall motion solve (JTJ) with the 
Gilbert form of damping^^^^^^^^ Thig ig gig _ 

nificant because, as we now discuss, Gilbert damping and 
Landau-Lifshitz damping produce quite different results 
for this problem when the same spin transfer torque is 
used. 

Consider a Neel wall where M lies entirely in the plane 
of a thin film when the current is zero. By definition, 
MxH = if we choose M(.t) as the equilibrium structure 
which minimizes the free energy .F[M]. The wall distorts 
if M x H 7^ for any reason. The theoretical literature 
cited above shows that, with damping omitted, the Neel 
wall moves undistorted at the speed v [see (O] when f3 = 



in ([2]). Gilbert damping brings this motion to a stop 
because Dq rotates M(x) out-of-plane until the torque 
from magnetostatic shape anisotropy cancels the spin- 
transfer torque. However, if the non-adiabatic term in 
([2]) is non-zero, steady wall motion occurs at speed f3v/a. 

Using this information, two recent experiments^^ 
used their observations of average domain wall velocities 
very near v to infer that (3 « a for permalloy nanowires. 
This is consistent with microscopic calculations (which 
include disorder-induced spin-flip scattering) that report 
(3 = a 29 or (3 rs oM for realistic band models of an itiner- 
ant ferromagnet. On the other hand, calculations for "s- 
d" models of ferromagnets with localized moments find 
little numerical relationship between (3 and a 30 - 31 . 

A rather different interpretation of the data follows 
from a discussion of current-driven domain wall motion 
in the s-d model offered by Barnes and Maekawa^ These 
authors argue that there should be no damping of the 
magnetization when a wall which corresponds to a min- 
imum of the free energy F[M] simply translates at con- 
stant speed. This is true of in because M x H = 
but it is not true of in ([5]) because M ^ when 
Nst ^ 0. From this point of view, the "correct" equa- 
tion of motion is 

M = - 7 M x H vd x M - AM x (M x H) , (6) 

because it reduces (for energy-minimizing walls) to 

M = -vd x M. (7) 

In the absence of extrinsic pinning, this argument iden- 
tifies the experimental observation of long-distance wall 
motion with a uniformly translating solution M(.t — vt) 
of with minimum energy. 

As we discuss below, it is possible to convert between 
descriptions with Landau-Lifshitz and Gilbert dampings 
by concurrently changing the value of the non-adiabatic 
spin-transfer torque. The Landau-Lifshitz description in 
Eq. [S] is equivalent to one with Gilbert damping with 
(3 = a. The goal of this paper is to argue that there are 
conceptual reasons to prefer the description with Landau- 
Lifshitz damping even when (3 ^ a. 

Section II presents micromagnetic simulations that 
confirm the discussion above and describes further de- 
tails. Then, the remainder of this paper provides three 
theoretical arguments which support the use Landau- 
Lifshitz damping for current-driven domain wall motion 
(in particular) and for other magnetization dynamics 
problems (in general). First, we reconcile our preference 
for Landau-Lifshitz damping with the explicit micro- 
scopic calculations of Gilbert da mping and non-adiabatic 
spin torque reported in Refs. l29H30H3ll . Second, we show 
that Gilbert damping can increase the magnetic free en- 
ergy in the presence of spin-transfer torques. Finally, 
we show that Landau-Lifshitz damping is uniquely se- 
lected for magnetization dynamics when the assumptions 
of non-equilibrium thermodynamics are valid. 
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FIG. 1: Position versus time for a transverse domain wall and 
several values of the applied current density computed with 
adiabatic spin torques (/3 = 0) and the two forms of damping 
in © and ©. 



II. MICROMAGNETICS 

Our analysis begins with a check on the robust- 
ness of the foregoing model predictions using full three- 
dimensional micromagnetic simulations of current-driven 
domain wall motion. 3 — We studied nanowires 12 nm thick 
and 100 nm wide with material parameters chosen to sim- 
ulate Ni8oFe2o- At zero current, this geometry and ma- 
terial system support in-plane magnetization with stable 
domain walls of transverse type^ Figure [T] shows the 
wall position as a function of time for a transverse do- 
main wall for several values of applied current density 
j. The curves labelled Gilbert (a = 0.02) show that 
wall motion comes quickly to a halt. Examination of the 
magnetization patterns confirms the torque cancellation 
mechanism outlined above. The curves labelled Landau- 
Lifshitz show that the wall moves uniformly with the ve- 
locity given by ([3]) which is independent of the damping 
parameter 

The sudden turn-on of the current and hence Oersted 
magnetic field at t = generates the small amplitude un- 
dulations of the curves in FigureQ]but otherwise has little 
effect on the dynamics. An initial state of a stable vortex 
wall in a 300 nm wide wire produces similar results, ex- 
cept that under the Gilbert formulation the vortex wall 
moves about twenty times farther before stopping as com- 
pared to the transverse wall in the 100 nm wire. We 
conclude from these simulations that the basic picture 
of domain wall dynamics gleaned from onc-dimcnsional 
models is correct. 

The magnetic free energy behaves differently in simu- 
lations depending on whether Landau-Lifshitz or Gilbert 
damping is used. Before the current is turned on, the do- 
main wall is in a configuration that is a local minimum in 
the energy. For Landau-Lifshitz damping, the energy re- 
mains largely constant near this minimum and is exactly 



2 



constant if the Oersted fields are ignored. For Gilbert 
damping, the energy increases when the current is turned 
on and the walls distort. For a transverse wall, the distor- 
tion is largely an out of plane tilting. Initially, the energy 
increases at a rate proportional to the damping param- 
eter (ignoring higher order corrections discussed in the 
next section) . The details of this behavior are somewhat 
obscured by the oscillations due to the Oersted magnetic 
field, but are quite apparent in simulations in which this 
field is omitted. As the wall tilts out of plane, the torque 
due to the magnetostatic field opposes the wall motion 
and the wall slows down. Eventually the torque balances 
the adiabatic spin transfer torque and the wall stops. 

In simulations using Gilbert damping, the change in 
magnetic free energy between the initial and final config- 
urations is independent of the damping parameter as it 
is determined by the balance between the magnetostatic 
torque and the adiabatic spin transfer torque. However, 
the amount of time before the wall stops and the distance 
the wall moves are inversely proportional to the damping 
parameter. The Gilbert damping torque is responsible 
for this increase in energy as can be seen from analyzing 
the directions of the other torques. Processional torques, 
like those due to the exchange and the magnetostatic 
interactions that are important in these simulations, by 
their nature are directed in constant energy directions 
and do not change the magnetic free energy. The adia- 
batic spin transfer torque is in a direction that translates 
the domain wall and does not change the energy in sys- 
tems where the energy does not depend on the position of 
the wall. Thus, in simulations of ideal domain wall mo- 
tion without Oersted fields, the Gilbert damping torque 
is the only torque that changes the energy. Throughout 
these simulations, the Gilbert damping torque is in a di- 
rection that increases rather than decreases the magnetic 
free energy. 



III. MAGNETIC DAMPING WITH 
SPIN-TRANSFER TORQUE 

When Nst = 0, it is well known that a few lines of 
algebra converts the equation of motion (TTJ with Gilbert 
damping into ([1]) with Landau-Lifshitz damping (and 
vice-versa) with suitable redefinitions of the precession 
constant 7 and the damping constants A and The 
same algebraic manipulations 39 show that ([6]) is math- 
ematically equivalent to a Gilbert-type equation with 
a = A/7: 



M = -7(1 + a 2 )M x H + aM x M 



v d r M - aM x d r M 



(8) 



To analyze (|8]), we first ignore spin-transfer (put v = 0) 
and note that this re-written Landau-Lifshitz equation 
differs from the conventional Gilbert equation only by 
an 0{a 2 ) renormalization of the gyromagnetic ratio. 



Consequently, first-principles derivations of any equa- 
tion of motion for the magnetization must be carried 
to second order in the putative damping parameter if 
one hopes to distinguish Landau-Lifshitz damping from 
Gilbert damping. This observation shows that papers 
that derive Gilbert damping2&2&2i>21i2&22 or Landau- 
Lifshitz dampin g 40 ! 41 from microscopic calculations car- 
ried out only to first order in a cannot be used to justify 
one form of damping over the other. 

Now restore the spin-transfer terms in ([5]) and note 
that the transformation to this equation from ^ auto- 
matically generates a non-adiabatic torque with = a. 
This transformation means that to lowest order in a and 
/3, an equation of motion with Gilbert damping and a 
non-adiabatic coefficient [3g is equivalent to an equa- 
tion of motion with Landau-Lifshitz damping with non- 
adiabatic coefficient = @g — a - This shows that 
equivalent equations of motion can be made using cither 
form of damping, albeit with rather different descriptions 
of current induced domain wall motion. Nevertheless, as 
we argue below, there are conceptual advantages to the 
Landau-Lifshitz form. 



IV. LANDAU-LIFSHITZ DAMPING UNIQUELY 
REDUCES MAGNETIC FREE ENERGY 

Landau-Lifshitz damping irreversibly reduces mag- 
netic free energy when spin-transfer torque is present. 
The same statement is not true for Gilbert damping. 
This can be seen from the situation described in Sec- 
tion II where Gilbert damping causes a minimum energy 
domain wall configuration to distort and tilt out of plane. 
Nothing prevents an increase in magnetic free energy for 
this open system, but it is clearly preferable if changes 
magnetic configurations that increase F[M] come from 
the effects of spin-transfer torque rather than from the ef- 
fects of a torque intended to model dissipative processes. 
This is an important reason to prefer in (|4|) to Dg 
in ||SJ| . This argument depends crucially on the fact that 
the adiabatic spin-transfer torque is not derivable from a 
free energy as we discuss henceforth. 

The field H in Eq. {J) is the (negative) gradient of 
the magnetic free energy. The component of this gra- 
dient in the direction that does not change the size of 
the magnetization is — M x [M x H]. Since this direc- 
tion is exactly that of the Landau-Lifshitz form of the 
damping, Eq. it follows that this form of the damp- 
ing always reduces this magnetic free energy. When 
the Gilbert form of the damping, Eq. is used in 
Eq. (TTJ, it is possible to rewrite the damping term as 
D G = -ajM x [M x H - (1/7)N ST ] + 0(a 2 ). Further, 
one can always write Nst = — 7M x Hst where Hst 
is an effective "spin transfer magnetic field" . However, 
unlike the field ^ H = -SF/SM in (JJ), there is no "spin 
transfer free energy" Fst which gives Hst as its gradient: 



Mo Hst 



SFi 



ST 



SM 



(not correct). 



(9) 
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FIG. 2: A one-dimensional Neel domain wall with magneti- 
zation M(x). 



If ([9]) were true, the lowest order (in a) Gilbert damping 
term -(17MX [M x (H+Hst)] would indeed always lower 
the sum F + Fst. Unfortunately, a clear and convincing 
demonstration of the non-conservative nature of the spin- 
transfer torque is not easy to find. Therefore, in what 
follows, we focus on the adiabatic contribution to ([2]) and 
show that a contradiction arises if © and its equivalent, 

dF ST = -MoH S t • dM, (10) 

are true. 

For this argument we consider a simpler model than 
that discussed in Section II. Figure [2] shows the magneti- 
zation M(x) for a one dimensional Neel wall in a system 
with uniaxial anisotropy along the x-direction. The do- 
main wall of width w is centered at x = and the plane 
of the magnetization is tilted out of the x-y plane by an 
angle A convenient parameterization of the in-plane 
rotation angle 8(x) is 

0{x) = tt/2 + sin _1 [tanh(x/u;)]. (11) 

Therefore, 

M = M[cos#(x),sin6>(x) cos0, sin0(x) sin^], (12) 

where cos#(x) = — tanh(x/w) and 

sin 6(x) = sech(x/w). (13) 

The magnetic free energy of this domain wall is indepen- 
dent of both its position and its orientation (angle 4>). 

For electron flow in the x-direction, |j2]) shows that the 
adiabatic piece of the spin-transfer torque lies entirely in 
the plane of the magnetization: 

oc 0'(x)(-sin0,cos0cos(?!>,cos0sin0). (14) 

This torque rotates the magnetization in a manner 
which produces uniform translation of the wall in the 
x-direction with no change in <f>. Since 

9'{x) = (l/w)sech(x/w), (15) 

comparison with (fl3|) shows that = outside the 

wall as expected. The magnetic free energy of the domain 
wall does not change as the wall is translated. 



Now, as indicated above ([9]), we are free to interpret 
the foregoing wall translation as resulting from local pre- 
cession of M(x) around an effective field Hst(x) directed 
perpendicular to the plane of the domain wall. Specifi- 
cally, 

Hst(x) oc 0'(x)(Q, — sin </>, cos </>). (16) 

However, if © and thus (flT)]) arc assumed to be correct, 
the magnitude and direction of Hgx imply that the puta- 
tive free energy Fst decreases when M(x) rotates rigidly 
around the x-axis in the direction of increasing </>.— On 
the other hand, the free energy must return to its original 
value when <fi rotates through 2ir. Since the gradient ([9]) 
can never increase the free energy, we are forced to con- 
clude that our assumption that Fst exists is incorrect. 

V. A LANGEVIN EQUATION FOR THE 
MAGNETIZATION 

Neglected work by Iwata^ treats magnetization dy- 
namics from the point of view of the thermodynamics 
of irreversible processes.— His non-perturbative calcu- 
lations uniquely generates the Landau-Lifshitz form of 
damping. In this section, we make equivalent assump- 
tions but go farther and derive an expression for the 
damping constant. Mori and co-workers did this using a 
projection operator method.— Our more accessible dis- 
cussion follows Rcif's derivation of a Langevin equation 
for Brownian motion.— 

We begin by taking the energy change in a unit volume 

dE = -fj, H a dM a , (17) 

where the repeated index a implies a sum over Carte- 
sian coordinates. It is crucial to note that the magnitude 
M = M is fixed so only rotations of M toward the effec- 
tive field H change the energy of the system. The interac- 
tion with the environment enters the equation of motion 
for the magnetization through a fluctuating torque N' a : 

^ = _ 7 (MxH) Q + <. (18) 

The torque N' is perpendicular to M since |M| = M. 

We consider the evolution of the magnetization over a 
time interval At which is much less than the precession 
period, but much greater than the characteristic time 
scale for the fluctuations r*. After this time interval, 
the statistical average of the change in magnetization 
AM a = M a (t + At) - M a {t) is 

rt+At 

AM a = - 7 (M x H) a (At) + / dt' < N' a {t') > . 

(19) 

The equilibrium Boltzmann weighting factor Wq gives 
< N' a {t') > = 0. However, < N' a {t') >^ when the 
magnetization is out of equilibrium. Indeed, this 
method derives the damping term precisely from the bias 
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built into the fluctuations due to the changes AE = 
—(XoHvAMv in the energy of the magnetic system. 

The Boltzmann weight used to calculate < N' a {t') > 
is W = W exp(-AE/(k B T)) where (assuming that H 
does not change much over the integration interval), 

AE(t') = -iM,H v tf) J g > dt" 

« -lM,H„{t)J N' v (t")dt". (20) 

Note that precession does not contribute to AE(t'). Only 
motions of the magnetization that change the energy of 
the magnetic subsystem produce bias in the torque fluc- 
tuations. Therefore, since W = W (l - AE/{k B T)) for 
small AE/ (k B T), the last term in (Til?|) now involves only 
an average over the equilibrium ensemble: 

AM a w - 7 (M x H) a (At) (21) 

+ dt ' f dt " < N ' a {t')K{t") > . 

We recall now that the torque fluctuations are corre- 
lated over a microscopic time r* that is much shorter 
than the small but macroscopic time-interval over which 
we integrate. Therefore, to the extent that memory ef- 
fects are negligible, we define the damping constant A (a 
type of fluctuation-dissipation result) from 

J dt" <N' a {t')N' v {t")> « XfaTM/tMdSi,, (22) 

for \t' - t\ > r* and with 6£ v = 5 alJ - M a M Vl which 
restricts the fluctuations to be transverse to the magne- 
tization, but otherwise uncorrelated. This approxima- 
tion reduces the last term in (|2"T|) to XMH± a At, where 
= — M x (M x H), is the piece of H which is perpen- 
dicular to M. Substituting (|22"|) into (f2~Tj) gives the final 
result in the form, 

— « - 7 (M x H) - AM x (M x H). (23) 
dt 

Equation (123|) is the Landau-Lifshitz equation for the sta- 
tistically averaged magnetization. It becomes a Langevin 
equation when we add a (now) unbiased random torque 
to the right hand side. 

The procedure outlined above generates higher order 
terms in A from the expansion of the thermal weighting to 



higher order in AE. The second order terms involve an 
equilibrium average of three powers of N' . These are zero 
for Gaussian fluctuations. The third order terms involve 
an average of four powers of N' , and are non-zero. They 
lead to a term proportional to X 2 H^H±, which we expect 
to be small and to modify only large-angle motions of the 
magnetization. 

VI. SUMMARY 

In this paper, we analyzed current-driven domain wall 
motion using both Gilbert-type and Landau-Lifshitz- 
type damping of the magnetization motion. Equiva- 
lent equations of motion can be written with either type 
of damping, but the implied description of the dynam- 
ics (and the relative importance of adiabatic and non- 
adiabatic effects) is very different in the two cases. 

With Landau-Lifshitz damping assumed, adiabatic 
spin transfer torque dominates and produces uniform 
translation of the wall. Non-adiabatic contributions to 
the spin transfer torque distort the wall, raise its mag- 
netic energy, and thus produce a magnetostatic torque 
which perturbs the wall velocity. Damping always acts to 
reduce the distortion back towards the original minimum- 
energy wall configuration. With Gilbert damping as- 
sumed, the damping torque itself distorts and thereby 
raises the magnetic energy of the moving wall. The 
distortion-induced magnetostatic torque stops domain 
wall motion altogether. Additional wall distortions pro- 
duced by non-adiabatic spin-transfer torque are needed 
to produce wall motion. 

In our view, Landau-Lifshitz damping is always prefer- 
able to Gilbert damping. When spin-transfer torque 
is present, this form of damping inexorably moves the 
magnetic free energy toward a local minimum. Gilbert 
damping does not. Even in the absence of spin-transfer 
torque, arguments based on irreversible thermodynam- 
ics show that the Landau-Lifshitz form of damping is 
uniquely selected for a macroscopic description. 45 Here, 
we proceeded equivalently and derived the Landau- 
Lifshitz equation of motion as the unique Langevin equa- 
tion for the statistical average of a fluctuating magneti- 
zation with fixed spin length. 
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Smith, G. Tatara, and Y. Tserkovnyak for useful discus- 
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